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Abstract 

Using the method of energy-level statistics, the 
localization properties of electrons moving in two 
dimensions in the presence of a perpendicular 
random magnetic field and additional random 
disorder potentials are investigated. For this 
model, extended states have recently been pro- 
posed to exist in the middle of the band. In 
contrast, from our calculations of the large- s be- 
havior of the nearest neighbor level spacing dis- 
tribution P(s) and from a finite size scaling anal- 
ysis we find only localized states in the suggested 
energy and disorder range. 

1 Introduction 

Many investigations, both analytical and numer- 
ical, have been undertaken to elucidate the lo- 
calization problem of non-interacting electrons 
in random magnetic fields (RMF) (see e.g., 
|, |, |, |]). The RMF model has been proposed 
for studying the quantum Hall effect with Lan- 
dau level filling v = 1/2 where the system of 
quasi-particles may be viewed ELS Si Fermi liquid 
like metal. In a mean field approximation these 
composite fermions |], [6|, |7]] are subject to a fluc- 
tuating local magnetic field with zero mean. 

However, up to now, no satisfactory agreement 
has been reached on the question, whether or not 
a metal-insulator transition exists in the thermo- 
dynamic limit. In recent numerical studies [[|, ||], 
current carrying states have been reported for a 
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RMF model including additional random disor- 
der potentials. The authors of Ref. || report 
that a metal-insulator transition has been de- 
tected in the band center at a critical disorder 
W c = 3.0 ± 0.2 by making use of the topological 
properties of the wave functions. 

The method of calculating the topological 
quantum number (Chern number) of the one- 
electron wave functions, which has proven to be 
a useful tool in the integer quantum Hall effect 
||, considerably restricts the achievable system 
size, L/a < 20. On the other hand, using the 
method of energy-level statistics [1C, |ll], 12, 13], 
system sizes of more than a factor of 10 larger 
are accessible. In this method, the knowledge 
of the eigenfunctions is not necessary, because 
the information about localized, critical or ex- 
tended states is gained from the spectral corre- 
lations. Therefore, we report here on an inde- 
pendent check of whether the extended states 
found in Ref. || can really be expected in the 
thermodynamic limit or if their result must be 
attributed to finite size effects. 

2 Model and Method 

A one-band tight-binding model is used to de- 
scribe non-interacting particles moving in two di- 
mensions in the presence of spatially fluctuating 
magnetic fields. 

H= £ Fexp(i&,)|fc><f|+£ e *l fe ><*l (1) 

<kj^l> k 

The \k) £11*6 £L complete set of lattice vectors con- 
nected with the sites of a square lattice with 



1 




Figure 1: The large-s behavior of the nearest 
neighbor level spacing distribution P(s) for dif- 
ferent system sizes L/a = 20 (□), 40 (o), 100 (o), 
and 200 (•). Disorder W/V = 2.75 and strength 
of the random flux ho = 0.5. 

linear size L and lattice constant a, and V is 
taken as the unit of energy. The sum runs 
over nearest neighbors only. The phases 4>kl 
along the bonds are chosen such that the cor- 
responding magnetic fluxes per plaquette, = 

4 > k,k+a x ~\~(Pk+a x ,k+a x +a y "'rtpk+ax+a.y ,k+a y 

are distributed according to — fiQir < < h^ir 
with probability density p(&k) = l/(2/io)- The 
strength of the flux was taken to be ho = 0.5. 

The diagonal disorder potentials are cho- 
sen at random from an interval [— W, W] with 
uniform distribution, p(ej.) = 1/(2W). We 
considered system sizes L/a in the range 20 
to 200 and disorder values W/V between 1.5 
and 3.5. The energy eigenvalues within a small 
range (|2£|/V < 0.3) in the middle of the tight- 
binding band are calculated using a Lanczos al- 
gorithm. For each pair of parameters (L, W) 
about 250 000 eigenvalues are accumulated tak- 
ing different realizations of the disorder poten- 
tials. We carefully checked that the results were 
independent of the width of the energy interval 
chosen. 

From the unfolded energy eigenvalues the 
nearest neighbor level spacing distribution P(s) 
is calculated, where s = l-Bj+i — Ei\/A and 
A is the mean level spacing. For infinite sys- 



tems it is known from random matrix theory 
(RMT) |L4], [TJ that in the metallic regime 
P (s) is close to the Wigner surmise, P(s) = 
Aas@ exp(— Bps 2 ), where f3 reflects the symme- 
try of the system: orthogonal (/? = 1), unitary 
(f3 = 2), and symplectic ((3 = 4). In the insulat- 
ing case the spacings of the uncorrelated eigen- 
values are well described by the Poisson proba- 
bility density distribution P(s) = exp(— s). 

In finite systems, when the correlation length 
exceeds the system size, one expects a hybrid 
of the Wigner and the Poisson form which even- 
tually tends to the limiting true metallic case, 
if W < W c , or to a complete insulating behav- 
ior, if W > W c , as the system size is increased, 
L — > oo. For disorder strength W close to the 
critical value W c various new scale independent 
distributions have been discovered recently in 
several systems [10, 11, [H| [u], 17] that exhibit 
a metal-insulator-transition or at least a critical 
point. 
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Figure 2: The calculated Jq(L, W) for various 
system sizes L and disorder strength W showing 
a behavior expected for localized states. 



3 Results and Discussion 

In the 2d random flux model with additional di- 
agonal disorder potentials, we find for small s, 
P(s) ~ s 2 , as expected for broken time rever- 
sal symmetry (unitary). Fig. [l| shows only the 
large-s behavior of P(s) for a disorder strength 
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Figure 3: The result of the finite size scal- 
ing analysis is a one-branch scaling function, 
Jo (£/£). All states are localized for disorder po- 
tential strength 2.25 < W/V < 3.5. 

W = 2.75 V, which is smaller than the proposed 
critical disorder W c = 3 V, and for system sizes 
L/a = 20, 40, 100, and 200. The RMT result 
for the Gaussian unitary ensemble (GUE) which 
holds in the metallic case and the Poisson result 
applicable for the uncorrelated eigenvalues of the 
localized states are also plotted. Increasing the 
system size shifts the curves away from the hy- 
brid P(s) towards the exp(— s) decay of the lo- 
calized states. Hence, even for W smaller than 
the suggested critical W c /V = 3, all the states 
still tend to localization in the limit L — ► oo. 

To be more quantitative, we have calculated 
the number Jq(L,W) = f^°Q(0,s)ds where 
Q(n,s) is the probability that a given energy 
interval s contains exactly n eigenvalues, and 
P(s) = d 2 Q(0,s)/ds 2 . In the metallic regime, a 
value Jq = 0.590 is known for infinite systems 
from RMT while for localized states Jg = 1. 
In Fig. |2] Jq(L,W) is shown for different sys- 
tem sizes and disorder strength. The values 
clearly increase with increasing L/a, at least for 
W/V > 2.25. This indicates localized states at 
the band center for this disorder range. A weak 
decrease of Jo (J, W) is still found for smaller dis- 
order strength, but no size dependence could be 
observed anymore, because in that case the local- 
ization length considerably exceeds our largest 



system size L/a = 200. 

It is possible to use Jq(L,W) also clS cl scal- 
ing variable, a procedure which was exploited re- 
cently [18 1 to determine the critical exponent of 
the correlation function at the metal-insulator- 
transition (MIT) in a 2d symplectic system. For 
the RMF model we have tested the validity 
of the single parameter hypothesis Jo (J, W) = 
/(L/£(W)) of a finite size scaling analysis which 
was performed for the data showing a size de- 
pendence (W/V > 2.25). The result is plotted 
in Fig. [| where only a single branch of the scal- 
ing function can be seen which excludes the exis- 
tence of a MIT for the disorder strengths inves- 
tigated. 




W/V 

Figure 4: Disorder dependence of the correlation 
function f (W) of the RFM (•) and of the 2d 
Anderson model (o) both normalized such that 
((W/V = 3.5) = 1. 

The disorder dependence of the scaling param- 
eter £(W) is shown in Fig. |I| in comparison with 
the result for a conventional 2d Anderson model 
without random magnetic fields. Since the pref- 
actor £o was n °t determined in our calculation, 
both curves were set to 1 at W/V = 3.5. The 
faster increase of £(W) in the RFM model indi- 
cates a weakening of the localization due to the 
random magnetic fields. The increase of the lo- 
calization length is, however, much larger than 
the factor of 2 expected from RMT [19] for the 
breaking of the time reversal symmetry. 
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In conclusion, the localization properties of 
a two-dimensional electronic model containing 
random magnetic fields and additional disorder 
potentials have been investigated numerically. 
Taking the same parameters as in Ref. || for the 
strength of the random flux, ho = 0.5, and the 
disorder potential W/V = 3, but system sizes 
up to L/a = 200, our data definitely show com- 
plete localization of the electronic states. Even 
for smaller disorder no sign of a delocalized phase 
could be detected. 
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